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Abstract 



To describe possible effects of medium, besides the hamiltonian introduced by 
Mikheev, Smirnov and Wolfenstein (MSW), in principle, there may exist new 
terms in the evolution equation for the neutrino oscillation due to the environ- 
ment. Considering a two-energy-level quantum system: ^1 — ^2, two-species 
neutrino system which is embedded in a matter environment, we solve the 
neutrino evolution equation precisely with the possible extra terms, especially, 
those cause alternations between pure and mixed states of the neutrino quan- 
tum system. We obtain some remarkable features for the neutrino oscillations, 
which arc different from that of MSW and accessible to be observed in solar 
neutrino observation and/or in the planned long-baseline neutrino experiments. 



PACS: 14.60Pq, 03.65Bz, 05.75.-b 



To solve the solar neutrino deficit problem Q , there have been many models which all assume 
neutrinos to be massive And it has been realized that the solar medium influence on the 
neutrino oscillation may be crucial to solving the problem. Due to the influence of the solar 
medium, the behavior of oscillation can be drastically modified and there can be a resonance in 
matter if sm9 m ~ 1, even though the mixing angle in vacuum 8 V is very small. This is the well- 
known Mikheev-Smirnov-Wolfenstein (MSW) mechanism H [|. Besides the solar matter effect, 
the Earth matter effect may also play a certain role in the neutrino oscillation problem, i.e. the 
recently reported atmospheric neutrino puzzle^) could be attributed to the medium effects either. 

The key point of the oscillation scenario may be illustrated by a two-level quantum system 
(y ei v^) or {yi,V2). The former are eigenstates of weak-interaction and the later are that of mass 
matrix respectively. In the scenario, the relevant quantities derived from quantum field theory are 
transformed into a Hamiltonian of quantum mechanics (QM). The species of neutrino produced 
by nuclear reactions is v e which can be decomposed with respect to the basis of v\ and v^. Thus, 
without the medium effects, its evolution is 

\v e {t) >= cos0„e -iElt |z/i > +sin6»„e~ iB2t |i/ 2 >, (1) 

where 9 V is the mixing angle in vacuum and E±, E% are energies of V\ and v%. As neutrinos are 
light, we have E{ w k + mf/2k (i = 1, 2), k = |k| and k is the three-momentum of the neutrinos. 

When MSW medium-effects are taken into account, the Hamiltonian should be modified from 
H to H' as 

H = k+(™ ° 2 \=>H' = k+( 4-GN e cos"6 v -GN e sm8 v co S v \ 
\ =J J \ -GN e sm6 v cos9 v ^-GN e sm 2 e v J 

where GN e is related to the electron density in a uniform matter environment. 

By solving the Schrodinger equation for the neutrino evolution 

one can see that a resonant structure of the mixing angle in medium 6 m is resulted in and the 
oscillation form becomes || 0] 

I < v»We(t) > | 2 = isin 2 (2^)(^) 2 [l -cos^], (4) 



where 

l m = l v (k)[l + (M^) 2 - 2co S (2^)(^)]- 1 /2, (5 ) 
to to 

and 

The mixing angle and period of the neutrino oscillation are modified accordingly. 

Recently, some authors have studied various extensions of quantum mechanics (QM) (or some- 
times called violation of QM), such as due to existence of 'micro black holes', a non-hermitian 
Hamiltonian is induced, thus a quantum system of pure state would turn into a mixed state, and 
CP, even CPT violation may be resulted^, [7]]. From a different motivation, Reznik studied a 
quantum system submerged in a 'proper surrounding environment' and he found that instead of 
a non-reversible evolution of the quantum system from a pure state to a mixed, there may be a 
periodic transition between pure and mixed states H. 

In fact, neutrino propagating in the dense medium of the Sun may be such a system. The 
oscillation between two species of the neutrinos, in principle, may be affected by the environment 
on its way to the surface of the Sun. First of all, we would like to consider the 'extension' by Reznik, 
and following him we write down a generalized Hamiltonian and find the solution of the evolution 
equation. One can note that it is indeed an extension of the Mikheev-Smirnov-Wolfcnstcin's result. 

In the following, we briefly recall Reznik's work and present our derivation of a modified ex- 
pression for | < v^\v e {t) > | 2 , as well as demonstrate how the modification practically affects the 
oscillation rate. 



Formulation (i) Unitary evolution between pure and mixed states for a quantum system em- 
bedded in a matter environment. 

The quantum system can be described by a density matrix p in general. If the system is in a 
pure state, one has Trp 2 = Trp = I, instead, when it is in a mixed state, Trp 2 < 1. In quantum 
mechanics, the Schrodinger equation is 

£p=[H,p]. (6) 
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In order to consider the environmental influence due to matter in very general cases, especially 
a possible unitary evolution of the system between pure and mixed states, Reznik proposed an 
operator p which relates to p, the reduced density matrix of the system, as 

P = pp\ (7) 

whereas pft p (in general ^ p) stands for the density matrix of the environment. At absence of an 
extra QM extension term, p still satisfies the same Liouville equation as (||). 
At the time t = 0, the system is supposed to stay in a pure state, so 

p(( = 0)E|^><^hp = pt. (8) 

Reznik shows that the equation (^|) can be generalized to 

ij t P = [H, p] + Lp + pR + gijKipK'j, (9) 

where L, R, , Ki and K'^ all are hermitian operators, gij are real parameters. 

It is noted that as long as the Kj in eq.(|9|) is not a unit matrix (do), i-e. it is not a trivial case, 
there cannot be a simple equation for the density matrix p. Thus people sometimes call the last 
term in eq.(|^) as an extension of QM, whereas being different in nature, a corresponding extra 
term is named as the violation term of QM in the case of ref. jFJ]. In the present case, rewriting 
eq.(||), we have 

ij f P =[H + L,p]+g ij (K i pK' j ft - h.c), (10) 

and in the last part of the equation, p does not appear as an independent variable directly at all 
(see below for more details). We can also easily prove that a resultant p = pp^ describes a mixed 
state as Tr(p 2 ) < 1 in the case. In fact, if there is no the last extension term, d/dt(Trp 2 ) = 
always holds, namely, the system remains in either pure or mixed state forever. It never transits 
from a pure state to a mixed or vice versa. Whereas when the 'extra' term exists, it is easy to 
prove d/dt(Trp 2 ) ^ 0, and it indicates that an alternation [J] between pure and mixed states occurs. 

(ii) Let us consider a system, that only two species of neutrinos are involved, thus they can be 
treated as a two-energy-level quantum system. 

1 Throughout the paper wc adopt the ward 'alternation' here instead of 'oscillation' adopted by [B, to avoid 
possible confution with the 'oscillation' for neutrinos. 



Following ref.|9| and considering the situation for neutrinos, we may write down the term H 
appearing in eq.(||) directly as 

H = k+( ^ ° 2 \ (11) 
V & ) 

in the basis of \v\ > and \v?, >. The trivial part A: is a unit matrix and can be omitted in our later 
calculations. After a transformation p — > p U with 



U = 



exp[-i J (R-H)dt'], 



according to ref.H the eq.(KJ) is turned to 



and 



i-^Hp + gtjKipK'j (12) 



H = H + L. 



The extra Hamiltonian L is hermitian which is induced by the environment, therefore can be 
absorbed into the original Hamiltonian H. 

Comparing eq.(|j|) and eq.(p"2|), H can be identified as H' and one can easily believe that the 
extra L is nothing new, but the Hamiltonian derived by Mikheev, Smirnov and Wolfenstein. Thus 

T _ ( -GN e cos 2 V —GN e sm9 v cos6 v \ . . 

\ -GN e sin 9 V cos 6 V -GN e sin 2 6 V ) ' ^ ' 

Therefore H = H + L is exactly the hamiltonian H' in eq. (||) . 

However, one can notice that there is one more extra term in eq.([T^), i.e. gijKipK'j which 
violates the hermiticity of p (but not p) and it would influence the behavior of the neutrino 
oscillation. For simplicity, later on we will write this term just as KpK' and the free parameters 
are absorbed in K and K' . 

(iii) H = H + L is a 2 x 2 hermitian matrix, so it can be decomposed as 

H = a a + aiCTi + a 2 (T2 + a 3 cr 3 . (14) 

The first term in this expression is trivial, so we will neglect it in later calculations, and furthermore 
from the explicit form of H given in eqs.(2) and (13), take ai = in later calculations. 
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Due to energy conservation, [K, H] = must be satisfied, a general form of K should be 

a(ai<7i + a 3 o- 3 ), 

where a is an arbitrary parameter. K' , being a hermitian 2x2 matrix, can be written as 

K' = Xx<Jx + A2CT2 + A3CT3, 

where a possible but trivial term, Ao<7o, has been dropped out. 

Just for illustration and convenience for later discussions, following ref.Q, we choose the sim- 
plest form for K' as A173 (i.e. Ai = A2 = and A3 = A). Thus the extension term becomes 



KpK' = (aiCTi + a 3 (73)/5Acr3, 



(15) 



where the arbitrary parameter a is absorbed into A. So far, we cannot derive this term from a 
well-established theory yet, so that we leave it aside as an open question. 

As usually adopted method |(], f?J [|, we write p in a four- vector form. In fact, this step is not 
trivial at all because of the QM extension term, (below we will give more discussions on this point). 
One can decompose p as 

P = PqO'o + P\<*\ + fao 2 + /5 3 cr 3 . 



Thus 



would become a form 



j t P = Hp + KpK 1 



i—Pa = (H%™ + 6H ab )p b , a,6 = 0,l,2,3, 



(16) 



where both H qm and SH are 4x4 matrices which correspond to H = H + L and KpK' respectively. 



Thus we 
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(17) 



It is easy to check that in the (4 x 4) form, Hu X 4,) = (H^ 1 + <5i/ a fc)4x4 still retains hermiticity. 

(iv) Solution of eq.(|l7j) can be obtained by diagonalizing the 4x4 matrix. Because it is a 
hermitian matrix, so it can be diagonalized via a unitary matrix V, as diag(Hu X 4\) = V H^^V^ . 



Thus 

Pa{t) = ]T V^ b V bc p c (t = 0)e-W (18) 

b,c 

where (5b are eigenvalues of Hu X 4j in cq.([l7|). Since at the moment t=0, the system is at a pure 
state, p(t = 0) = p(t = 0). Without losing generality we assume a pure v e at the initial state, thus 

p(t = 0)=p(t = 0) = \v e Xv e \=( C °f 6 " -^Yr 6v )- (19) 
ry ' ry V - sm ^ cos ^ sm^„ / v ; 

Converting the initial condition in 2 x 2 matrix form into a four-vector form and substituting those 

p c (t — 0) into the solution fll8|), we obtain p a (t) and p e (^) = Pe(t)f>t(t) in (2 x 2) matrix form. 

In fact, a complete form of KpK' = a(ai(Ti + a3£73)p(Ai<7i + A2O2 + A3O3) which replaces the 
simple one in eq.(p"5|), does not change the whole physical picture. Our full derivation shows that 
to obtain the formulation with all the A; (i=l,2,3) existing, one only needs to re-write A of eq-dl^) 
in a complicated combination of A^ whereas the general form remains unchanged. So a solution 
derived from the simple expression does not lose generality. Here we omit the details to save space. 

(v) The oscillation rate is defined as 

P„.^„„(t) = I < u^\u e (t) > I 2 = Tr(p+p e (i)). (20) 

Thus we have the final expression of P as 

P(0 = -sin 2 20„(-p) 2 [l-cos( — )cos(— - ) + cos 20, sin( — )sin(-— )], (21) 

where Z is the distance from the production site of v e to the detection point. Note that our A is 
dimcnsionless and it is slightly different from Reznik's notation where his A has a dimension of 
energy. It is a modified expression of the neutrino oscillation in medium and simply an extension 
from that given by Mikheev, Smirnov and Wolfenstein. 

Note that we have checked that the extension of QM which we are considering, is of unitarity. 
Because of conservation of possibilities the constraint 

is satisfied exactly. 

As A = 0, which implies the extra term KpK 1 disappears, the expression reduces back to the 
form given in ref.(3) and ||. For non-zero A values, the oscillation is not simply harmonic, but in 
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a modified mode. In general the oscillation period is shortened, so for the solar neutrino case the 
result favors a larger reduction rate of the electron neutrinos produced in the Sun. 

(vi) To compare various possibilities, let us consider another type of environmental influence 
on the neutrino oscillation due to the 'micro black holes everywhere' although the authors of ref.^J 
originally proposed it for the quantum system K° — > K° mixing. According to the mechanism 
given in ref. S , the oscillation rate may be obtained: 

P{t) = | < f>e(*) > I 2 = ^sin 2 2^(l -e^^cos^), (22) 

where an approximation 

^ > \a\, \p\, | 7 | (or a ~ (5 > |/3|) 

is required. 

That is a modified vacuum oscillation and the damping factor indicates an energy loss at the 
evolution process. It is understood that as the neutrinos propagate in the physical vacuum full 
with micro black holes, their evolution behavior are affected by the black holes. Simultaneously the 
quantum system turns from a pure state into a mixed state and this trend is non-reversible since 
the black holes never release anything out (here we do not refer to the quantum tunneling effect 
of black holes) . The new behaviors of the system evolution will be discussed in our coming work []l3[ . 

Discussions and implications (i) It is definite that the solar environment can significantly 
influence the neutrino oscillation. Mikhcev, Smirnov and Wolfenstein's work accounted such effect 
by adding an extra hamiltonian to the system. Here in a more general framework, as Reznik 
suggested, there is a natural way to include the influence from environment, especially the system 
may be oscillating between pure and mixed states. Namely, the environment interacts with the 
quantum system in an unknown and complicated way. Such an interaction is not only depicted 
by a hermitian Hamiltonian, but also by an additional term, which violates the hermiticity of p. 
It makes the system evolve between pure and mixed states. The MSW's Hamiltonian only results 
in a hermitian part of p. It is also noted that the interaction of the neutrino system with the 
magnetic field of the Sun p~0) can also induces an external environmental effect but it may also 
be absorbed into L. Only the extension part can cause the system to alternate between pure and 
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mixed states j9|. 



Here the situation is different from that of the 'micro black holes'. There are interactions 
between the neutrino system and the environment, although we cannot theoretically derive them 
precisely at present, In general the energy of the whole system is conserved, even though it evolves 
from a pure state into a mixed state and vise versa. From the analytic expression we also 
observe that the oscillation behavior of the neutrino system with the QM extension term is quite 
different from that without it. In the regular QM region where KpK' is absent (i.e. A = 0), the 
oscillation rate | < v^\v e {t) > | 2 can vary between and sin 2 (2# m ) (or sin 2 (20 lI ) for the vacuum 
case) H ||], on the contrary, due to the presence of the QM extension term, | < Ufj,\u e (t) > | 2 can 
never be zero except at the moment t=0. 

For the solar neutrino problem since the initial v e is produced all over the Sun and the distri- 
bution of the chemical elements in the Sun varies, one must integrate the transition probability 
over the whole production region and energy spectrum of neutrinos. It is 



where the origin is chosen at the center of the Sun, and the time, which the neutrino takes for 
traveling from I to the surface of the Sun, is t = l/v ~ l/c and v is the speed of the neutrino, which 
is close to that of light c, as long as neutrino masses are supposed to be very tiny. The factor 1/6 
accounts for an average probability in the propagation direction to the Earth. If the distribution 
of the neutrino production rj(L — l,E) does not vary drastically, it can be approximately treated 
as a constant, the integration over sites can be carried out analytically. 

To solve the solar neutrino problem with MSW solution as suggested by many literatures, 
A = m 2 — m 2 , is comparatively large. We have L ^> i TO |0|, i.e. the solar size is much longer than 
the oscillation length. In the case the extra contribution in eq.(21) would disappear as well as the 
oscillation term in cq.(4) which is proportional to |% cos(2irL/l m ), unless A is close to unity. If 
A ~ 1, the modified expression (21) can be averaged as 



Note that, as long as A ^ 1 and l m <C L, the extension term KpK' does not change the physics 
obtained in the regular framework without such a term, whereas A — > 1, an extra factor cos 2 9 V 





(23) 
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would be attached to the original expression that is observable. If 6 V is small, this effect can also 
be neglected. 

In fact, the medium is not uniform, i.e. rj(L — I, E) is not a constant over I and A may be 
environment dependent, there may be some effects which can influence the data fitting of 6 m , but 
to take into accounts of all the detail, one needs better understanding of the solar model and the 
QM extension effect. 

It is also interesting to note that if A = m 2 — m\ is very small, for instance, which is close 
to the expected value for the vacuum oscillation solution, i.e. the oscillation in medium cannot 
complete in one cycle in the solar neutrino case, the effects of the extension term, no matter from 
the micro black holes or the dense environment, increases the transition probability from v e to v^. 

(ii) The interesting QM extension effect is possibly observable in other neutrino oscillation 
experiments. The most likely possibility is in the planned long baseline accelerator experiments, 
KEK-SuperKamiokande(250Km), CERN-Gran Sasso(730Km) and Fermilab-Soudan2(730Km)Jl|. 
The average energies of fi neutrino beams are approximately 1 GeV, 6 GeV and 10 GeV respectively. 
In all these long baseline experiments neutrino oscillation is expected to be seen if the neutrino 
mass square difference is not much smaller than 10~ 3 eV 2 and the mixing angle is not too small. 
With the new QM extension effect being taken into account, we need to consider two cases. If A 
is not much smaller than one, then the new effect may modify neutrino oscillation significantly. 
And if A is much larger than one, then the oscillation is observable for even smaller neutrino mass 
square difference 10 _3 A _1 eV 2 . If the distance of the experiment is chosen as 10 4 Km, that is 
about the diameter of the Earth, the energy of the neutrino beam is about lGeV, for neutrino 
mass square difference as small as 10~ 5 eV 2 , the oscillation is expected to be seen. So long as A is 
not too smaller than one, then the QM extension effect is also to be observed for the mass squared 
difference not too smaller than 10~ 5 eV 2 . 

(iii) Ellis et al. || and Reznik || introduced two different mechanisms which both can induce 
transition (or oscillation) of quantum systems from a pure state to a mixed one. Ellis et al. added 
the extra 5H^ x4 ^ of (4 x 4) form in the evolution equation for p v of the four- vector form and 
this new term is not hermitian. When one tries to turn the equation back to a (2 x 2) form, it is 



9 



found that there is no way to make a Schrodinger-type equation for the (2 x 2) p— matrix such as 
i(d/dt)p = H(2x2)P- Reznik's formulation, instead, gives an equation for p of 2 x 2 form where he 
introduced the extension term KpK' and the equation also cannot be converted into an equation 
for p, even in the four- vector form, but there Huxi) is a hermitian matrix. Generally the resultant 
(2 x 2) p is not hermitian, while p is. One can notice that the differences of Ellis et al.'s mechanism 
from Reznik's explicitly. 

Most generally, if p (or p) is in a four vector form, the generalized Schrodinger-type equation 
reads 

i^_P V = H {ixi) p v ', (or for p), 
it can be converted into a (2 x 2) matrix equation where p (or p) is a 2 x 2 matrix as 

i^p = Lp + pR + LpK + Ap, (or for p), (24) 

where L, R, K, A are operators and can be decomposed as 

O = O o- + O l a l (O = L, R, K, A; i= 1, 2, 3) , 

the coefficients O a (a — 0, 1, 2, 3) would be complex as the Lp(p)K term exists. Writing out iJ(4 X 4) , 
the real part of O a correspond to the hermitian components of H(axa) while the imaginary parts 
to the anti-hermitian ones. 

Because of the existence of the anti-hermitian components, the energy of the system is not 
conserved as in Ellis et al.'s Hu X 4\. In our formulation, the i?(4 X 4) is hermitian, so the energy of 
the system is conserved. Both the extra term KpK' and KpK' demand (d/dt)Tr(p 2 ) ^ 0, i.e. a 
system evolves from a pure state into a mixed state. 

(iv) We investigate the general properties of a quantum system which evolves from a pure state 
into a mixed one due to either the micro black holes in the background or certain environment in 
the case of neutrino oscillation. With the generalized Schrodinger equation, we have derived an 
expression for p which is non-hermitian and the density matrix p = ppK In this result, p ^ p\ 
naturally, one can decide that p 1 p describes the environment In the case, (d/dt)Trp 2 ^ 0, thus 
the system evolves between pure and mixed states. 
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Thus we have derived the v e — > transition rate and find that the newly obtained expression 
is a modification of the famous Mikheev-Smirnov-Wolfenstein's formulation. 

In the modified expression, sin 2 (20 t ,)(j a -) 2 = sin 2 (26> m ) keeps unchanged as in the MSW the- 
ory, but the oscillation form (1 — cos(y^)) is modified into a more complicated one as (1 — 
cos(2^) cos(^) + cos(20„) sin(^) sin(^)). It obviously favors v e -> v^. 

(v) The new extra term(s) either from micro black holes or from a matter environment should 
be derived from certain theories in principle, whereas so far no one has a reasonable way to derive 
it (them) out theoretically. However as pointed out in the paper, the parameters, such as A for 
instance, can be phenomenologically determined by fitting data if they fall into a suitable region. 
The environmental influence may be observable in the solar neutrino and/or the planned long- 
baseline experiments through neutrino oscillations. Therefore we conclude that the environmental 
influence may be experimentally accessible and certain constraints on the 'extra' terms may be set 
from the solar neutrino observation and the planned long-baseline neutrino oscillation experiments 
in the foreseeable future. 
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